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Symmetry-breaking phase transitions may leave be- 
hind topological defects [Q] with a density dependent on 
the quench rate We investigate the dynamics of such 
quenches for the one-dimensional, Landau-Ginzburg case 
and show that the density of kinks, n, scales differently 
with the quench timescale, tq, depending on whether the 
dynamics in the vicinity of the critical point is overdamped 
(n oc Tq^^"^) or underdamped (n oc Tq^^^). Either of these 
cases may be relevant to the early Universe, and we de- 
rive bounds on the initial density of topological defects in 
cosmological phase transitions. 
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Dynamics of symmetry breaking phase transforma- 
tions is of interest in cosmology, condensed matter and 
high energy physics. Recent surge of interdisciplinary 
interest has been fueled by the experiments involving 
creation of topological defects during rapid second or- 
der phase transitions in superfluids and other sys- 
tems iQ , in a setting reminiscent of similar cosmological 
processes [0,D- 

In the cosmological context, topological defects such 
as cosmic strings may have played a role in seeding 
structure formation |^,^. In high energy physics, ac- 
celerator experiments may allow one to probe restora- 
tion of some of the symmetries (e.g. chiral), which 
were originally broken early in the history of the Uni- 
verse. The signature of whether such restoration has 
occurred will come from the fluctuation of decay prod- 
ucts ||l^,|ll|, determined by the relevant dynamics of 
the order parameter during a quench. In superflu- 
ids, the very same process controls the production of 
topological defects. The interest in experimental explo- 
ration of the critical dynamics is therefore well-justified 
by its wide-ranging applications, which may also come 
to include in the near future the creation of vortices in 
Bose-Einstein condensates ]T^ . 

These processes span many orders of magnitude in 
spatial and energetic scales. Yet, a large class of them 
is well-approximated by the Landau-Ginzburg theory. 
Therefore, the dynamics of the order parameter (p is 
governed by an equation of the form : 

(p + r)(p- c^V^ip +[P(p^ - e{t) (p]/2 = '!9(i, x) . (1) 

Above, rj characterizes viscosity, while c, (3 and m are 
constant coefficients, and e{t) is the time-dependent 
relative temperature, assumed to vary with time as 



e ^ t/rq , where tq the is quench timescale. The term 
t9(i, x) is noise characterized by its spatial and temporal 
correlations, as well as by its amplitude 9. We assume: 
(i}{x,t),^{x',t')) =2T]eS{x'~x)S{t'-t). Equation (|l]) 
can be expressed in "natural" units t — > t/m, x 



x c/m, r] rim, ip 
which leads to: 



ipm/^/p and 



(fi + 7](p - 
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Thus, in the vicinity of the second-order phase transi- 
tion, an enormous range of "bare" parameters can be 
reduced to two: the "renormalized" damping rate rj 
and the noise temperature 9. The quench adds a de- 
pendence of the consequences of critical dynamics on 



the quench rate e 



'Q 



The aim of our study is 



to investigate the dependence of the critical dynamics 
on the value of 77 and tq, under the assumption that 
9 is sufficiently small, so the probability of thermally 
activated symmetry restoration is negligible after the 
quench. In this paper we take 9 — 0.01. We focus on 
creation of kinks — one-dimensional topological defects 
— in the course of rapid quenches. 

Evolution generated by Eq. (^ is overdamped when 
rjif > Lp. In this regime, the relaxation time ~ 
\'p/'p\ scales with the relative temperature e as: t,^ ~ 



Vrl\e\ 



^ ~ Tyrg |t| ^ , in the units of Eq. ([|) with 



To = TO ^ the dynamical timescale. In accord with |^, 
one expects the initial size ^ of the pieces of the new 
broken symmetry phase to be set at the time t, when 
the time to (from) the phase transition is comparable 
to the relaxation timescale and the freeze-out of 
the field evolution occurs; that is, its state cannot keep 
up with the change of the thermodynamic parameters 
as a result of critical slowing down. This freeze-out 
condition, T^{t^) = t^p yields 



U ~ To (tITq) 
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The correlation length which sets the stage for the 
defect formation Q|, is then: 
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where = c/m characterizes the low temperature (e 
0) healing length. 
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To test these arguments , we have recently carried 
out a numerical study of defect formation , showing 
that the density of kinks formed in a quench indeed 
varies, in the overdamped regime, as: 
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with / ~ 8. Similar study was independently carried 
out by Lythe who has estimated / — 27r(ln0)^/'* 
when 9-^1. 

Our purpose here is to extend these studies from the 
regime where damping dominates (which is most rel- 
evant in condensed matter applications) to the range 
where the evolution is underdamped (as may be the 
case in cosmology) . For details of the numerical tech- 
nique, see Ref. 0. 

In the underdamped case, (/3 will dominate, and 
the order parameter reacts to the quench-induced 
changes in the effective potential on the timescale 
~ \(p/(p\^/^. Thus, ~ To|e|-i/^ The freeze- 
out condition, T^{t^) = t(p yields in this underdamped 
regime: 



(7) 
(8) 



Consequently, the scaling of the characteristic correla- 
tion length with the quench rate Tq^ is expected to 
change to 



1/3 



(9) 



Furthermore, the density of the number of kinks is 
given in this case by: 
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(10) 



although / may be now different. 

We can therefore draw two related conclusions: (i) 
In the overdamped regime, the density of kinks should 
scale with r/^^^, and should become viscosity indepen- 
dent in the underdamped case, (ii) Power-law depen- 
dence of the density of kinks with the quench timescale 
should change from cx Tq in the overdamped case 

— 1/3 

to (X Tq in the underdamped case. The overdamped 
scalings should apply when the evolution is dominated 
by the first derivative {rjip > (f, i.e. > 1/t|) at 

the instant when topological defects "freeze-out." This 
will happen for: je^j > |e<p| , or — using Eqs. (^) and 
(^ — when: 



iVTof > (to/tq). 



(11) 



We identify kinks as zeros of the order parameter. 
This can be justified only well after the phase transi- 
tion, when if has locally settled into the broken symme- 
try state. Kinks annihilate, and their number slowly 



decreases with time. Previously, in the overdamped 
regime, we were able to confirm the predicted de- 
pendence of the initial number of kinks on tq from the 
numerical data by using a fairly straightforward proce- 
dure of simply counting zeros at a fixed value of I/tq 
[0. The nature of that dependence did not change 
dramatically even when the counting of kinks was tak- 
ing at a constant value of t (although a change on the 
slope as well as evidence of the saturation in the num- 
ber of kinks for small tq were noted). But the nature of 
critical dynamics and especially the annihilation rate 
depend on 77, which we shall vary by several orders 
of magnitude. To compare "initial densities" of kinks 
now, we therefore need a more objective procedures in- 
dependent of the time at which the kinks are counted. 
We have done this by using whole runs of kink densities 
(such as the ones shown in Fig. ^ to model annihilation 
either as a power-law N ~ No{t / tq)~^ , or as an expo- 
nential, N ~ iVoCxp {—at/TQ), with N the number of 
zeros of the order parameter. The actual dependences 
are usually sufficiently similar to a straight line that 
both of these procedures yield reasonable initial kink 
numbers No, and at least in some cases an almost sur- 
prisingly similar dependences on the parameters (see 
Fig. I). 

The dependence of the initial number of kinks No 
on the damping coefficient 77 for three different quench 
rates (tq = 128, 256, and 512) is shown in Fig. |^. In 
the regime of large viscosities, critical dynamics is over- 
damped, in accord with Eq. (pi]), and leads to the 
power-law dependence cx 77^/*, Eq. (||). The spacing be- 
tween the three lines is also roughly consistent with the 
one anticipated from that equation. As the damping 
rate decreases below the value estimated from Eq. ( |ll| ) , 
namely r] ^ 0.1, the number of kinks becomes essen- 
tially independent of rj. Moreover, the spacing between 
the (now approximately horizontal) lines of constant tq 

— 1/3 

is consistent with the underdamped case No cx Tq , 
Eq. (|0l). 

We should note, however, that for 77 < 10^^, the 
number of kinks are small and the annihilation is more 
efficient. Consequently, our results in this range are 
less reliable. In a sense. Fig. ^ indicates a more sys- 
tematic (and probably more consistent with the the- 
oretical expectations) trend with the damping rate 77 
than our estimates of the number of kinks plotted in 
Fig. ||. Most of the scatter in Fig. ^ stems from our 
inability to model annihilation rate in a consistent fash- 
ion over a broad range of parameters, rather than from 
the "raw" data shown in Fig. |l|. 

These conclusions concerning the transition from 
overdamped to underdamped behavior are strength- 
ened by comparing families of simulations correspond- 
ing to the same 77 but for varying tq (see Fig. ^ and 
As before, we consider two methods for obtaining the 
initial number of kinks Ng from the time-dependent 
data. Fig. ||. In the range of long quench timescales 
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they produce similar (but not identical) power-laws. 

According to condition (11), the 77 = 5 and 1 cases 
in Fig. ^ are, for the values of tq under consider- 
ation (2 < tq < 4098), entirely within the over- 
damped regime. For these two cases, we find power- 
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consistent with Eq. 



On the other 



laws . 
hand, for the rj — 1/5 case in accord with Eq. (|l|), 
a transition between overdamped and underdamped 
regimes should occur at tq 125. We find (see Fig. ^ 
indication of a change in the power-law dependence, 
from No oc Tq^^^ to oc Tq^^^ as rj decreases. 

One obvious case of breakdown of power-laws, 
Eqs. (^) and (|o|), occurs when the condition e <C 1 
of applicability of the theory of Ref. |^ is not satis- 
fied [|^. In that case, the predicted initial separation 
of kinks would be comparable to (or even smaller than) 
the zero-temperature healing length fo. This would of 
course result in a rapid initial annihilation, so that the 
density of defects would be set by the annihilation pro- 
cess rather than by the critical dynamics. We have seen 
evidence for this behavior for sufficiently small tq in 
Ref. 0. 

Equations (^) and (^^ for the initial density of topo- 
logical defects can be used in the cosmological setting. 
Phase transitions are likely to occur in the radiation 
dominated era, when the temperature T of the plasma 
and the Hubble time tn since the Big Bang are tied 
with the equation =^ constant. This immediately 

yields quench timescales tq = 21^ = H^^ , where H 
is the Hubble parameter. 

Damping rate is the other important parameter set 
by cosmology. In the radiation-dominated epoch -q = 
3 -I- 7 , where 3 iJ is the effective viscosity caused by 
the Hubble expansion, while 7 is damping due to the 
coupling with the other degrees of freedom. Early on, 
the "Hubble viscosity" may even dominate. 

The nature of the critical dynamics in the immedi- 
ate vicinity of the phase transition is decided by the 
inequality (pT|), which now reads 
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In the overdamped case 



e^ToH{3 + 7/H) 
This in turn leads to 



1/2 



(12) 

(13) 
(14) 

(15) 



where we have set c — 1, so — m^^. Thus the 
density of topological defects is principally set by the 
geometric average of the characteristic length scale of 
the order parameter (^o) and the size of the horizon 
(H^^). For small ^/H, damping is dominated by Hub- 
ble expansion. In that regime, H/m > 1 (or H^^ < ^o) 



would be required for the critical dynamics to be over- 
damped. This would lead to ^ ^ H^^, which in effect 
implies that in this case the density of defects is set by 
the size of the horizon. 

When the critical dynamics is underdamped; 



£~T„(roi/)-i/^ 



This immediately leads to 

i 0, ^o {H-^i^oY" 



(16) 
(17) 



(18) 



The characteristic distance ^ is then bigger than the 
healing length by the third root of the size of the 
horizon at the time of the transition measured in the 
units of ^o- 

We conclude by noting that the dependence of the 
number of kinks on the viscosity parameter corrobo- 
rates anticipated existence of the two regimes in the 
critical dynamics, each with a distinct scaling of the 
relevant characteristic timescale with the relative tem- 
perature e. Overdamped regime produces kinks with 

separations ^ oc {tq/t])^^^, while in the underdamped 
1/3 

case ^ oc Tq and is independent of viscosity. The 
subsequent annihilation rate of the kinks strongly de- 
pends on viscosity, and is much more rapid in the un- 
derdamped case. The borderline between the two is 
consistent with the considerations of |p|,|l7|. 
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FIG. 1. Number of kinks A*' as a function of time 
for simulations with different viscosity parameters r], but 
with the same quench timescale tq — 256. Two mod- 
els for the annihilation rate are shown, (a) exponential 
N = No exp {-at/rq) and (b) power-law = No{t/TQ)~^. 
Note the increase of the annihilation rates for small rj. 





FIG. 2. Initial number of kinks No as a function of 
the damping rate r; for a fixed quench rate timescales (top 
to bottom TQ — 128, 256 and 512). Both exponential (a) 
and power-law (b) model results (see Fig. |l|) are shown 
(and are essentially identical). For rj > 0.1, No oc r;®, 
where No is obtained from the fittings in Fig |^. Best 
fits yield g = (0.27 ± 0.035, 0.25 ± 0.029, 0.27 ± 0.011) for 
TQ = (128, 256, 512), respectively. 




FIG. 3. Number of kinks A as a function of time for a 
fixed viscosity rj = 1 but different quench timescales tq. As 
with Fig ^, (a) corresponds to an exponential fit and (b) to 
a power-law fit. Saturation is apparent in the case (a) for 
short quench timescales (and large initial kink densities). 
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FIG. 4. Dependence of the initial number of kinks 
Nc on the quench timescale tq for values of the damp- 
ing rate 77 = 5 (top), 1 (middle) and 1/5 (bottom). 
Case (a) is obtained from an exponential fitting to the 
decay of number of kinks and (b) from a power-law 
fit. Fittings to No oc Tg* yield (from top to bottom) 
g = (0.23 ± 0.010, 0.26 ± 0.011, 0.33 ± 0.011) in case (a) 
and g = (0.28 ± 0.010, 0.30 ± 0.011, 0.36 ± 0.010) in case 
(b). 
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